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L^i ABSTRACT 
> 

t^j- ■ We present the chemical compositions for eight bright giants in the globular 

cluster NGC 1851. Our analysis reveals large star-to-star abundance variations 
and correlations of the light elements 0, Na, and Al, a feature found in every 
well studied globular cluster. However, NGC 1851 also exhibits large star-to- 
^ ■ star abundance variations of the s-process elements Zr and La. These s-process 

elements are correlated with Al, and anticorrelated with 0. Furthermore, the Zr 
and La abundances appear to cluster around two distinct values. A recent study 
^ ■ revealed a double subgiant branch in NGC 1851. Our data reinforce the notion 

that there are two stellar populations in NGC 1851 and indicate that this cluster 
has experienced a complicated formation history with similarities to u Centauri. 



Subject headings: Galaxy: Globular Clusters: Individual: NGC 1851, Galaxy: 
Stars: Abundances 



1 Based on observations made with ESO Telescopes at the Paranal Observatories under programme ID 
70.B-0361(A) 
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Introduction 



Recent studies have revealed that the Galactic globular c lusters u Ce ntauri (IBedin et a 



20041) . NGC 2808 flD'Antona et al.ll2005l : IPiotto et al.l 120071 ). and M 54 (ISiegel et al.l 120071 ) 



contain multiple stellar populations. The evidence consists of HST color-magnitude dia- 
grams (CMDs) in which multiple main sequences have been identified. In all cases, the 
most plausible explanation is that these globular clusters contain populations with distinct 
compositions and/or ages, in contrast to the classical concept that globular clusters are 
monometallic, coeval, simple stellar populations. The most recent addition to this intriguing 
collection is NGC 1851, a globular clust er with a bimodal horizontal branch (HB), whose 
CMD displays a double subgiant branch (IMilone et al.l 120071 ). 



Both u Centauri and M 54 possess a large metallicity spread and are the two most mas- 
sive Galactic globu lar clusters, a; Ce ntauri is widely regarded as the nucleus of an accreted 
dwarf galaxy (e.g., ISmith et al.l 120001 ) and M 54 is suspected to be the nucleus of the Sagit- 



tarius dwarf spheroidal galaxy (e.g., lLayden fc Sarajedinill2000l ). However, not all clusters 



with multiple main sequences show a metallicity spread and although there i s speculatio n 
that all globular clusters may have been the nuclei of ancient dwarf galaxies (IBekkil 120061 ) . 
currently there is little direct observational evidence. For NGC 2808, the multiple main 
sequences and complex HB morphology can be attributed to a larg e He abundance variation 



lp 

flD'Antona k Calol [20041 : iD'Antona et al.ll2005l : IPiotto et al.ll2007D . but the Fe abundances 
and ages show little, if any, variation among the populations. Although ther e are large star- 
to-star abundance variations of O, Na, and Al in NGC 2808 ( Carrettal 2006), such pa tterns 



are ubiquitous in Galactic globular clusters (e.g.. lKraftlll994l and lGratton et al.ll2004l ). 



The situation regarding the chemical compositions for stars in the globular cluster NGC 
1851 is less cl ear due to the lack of data. The sole spectroscopic abundance analysis was 
conducted by iHesser et al.l (119821 ). Low resolu tion spectra r eveale d extremely strong CN 
bands for three out of eight bright red giants. IHesser et al.l (1l982f ) also found that the Sr 
and Ba lines may be enhanced in CN-strong stars relative to other giants in this cluster. 
Such abundance patterns, if confirmed, would bear a striking resemblance to uj Centauri. 
A study of the stellar chemical compositions in NGC 1851 is of great interest in order to 
provide constraints upon the interpretation of the double subgiant branch and to explore 
the chemical evolution of this cluster. 
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2. Observations and abundance analysis 



The spectra for eight bright giants in NG C 1851 were retriev ed from the ESO archive. 
The observations were obtained using UVES (IDekker et al.ll2000l ). Each star was observed 
for 900 seconds using the 0.70" slit which provided a spectral resolution of R ~ 55,000 per 
4 pixel resolution element. The data were reduced using standard procedures in IRAF. The 
wavelength coverage ranged from 5950A to 9750A with a typical signal-to-noise ratio of 60 
per pixel near 6500A. All stars are radial velocity m embers with y r rang ing from +316 km/s 
to +332 km/s. One star, 333, was also observed by lHesser et al.l (119821 ) . 



The stellar param eters were determined u sing the same tools and techniques described in 



Grundahl et al.l (120021 ) and lYong et al.l (120061 ) . While the agreement between the photomet- 



ric and spectroscopic stellar parameters was excellent, we adopted the spectros copic param- 
eters. (Photometric e stimates were based on data from iGrundahl et al.l (Il999l ) and 2MASS 
(ISkrutskie et al.ll2006l ). The results and conclusions would not change had we adopted pho- 
tometric stellar parameters.) 

Abundances were measured from an equivalent width analysis. For O, La, and Eu, 
synthetic spectra were generated to account for blends, hyperfme structure, and/or isotopic 



splitt i ng. We used th e LTE line analysis and spectrum synthesis program MOOG (Sneden 



9731) . iKuruca (119931 ) model atmospheres, and the same line lists employed by lYong et al 



( 20051 ). While many lines of interest lie below 6000A and were not observed, there were a 



sufficient number of lines from which the abundances of O, Na, Mg, Al, Si, Ca, Sc, Ti, V, 
Mn, Fe, Co, Ni, Zr, La, and Eu could be measured. We focus upon the most fascinating 
results which come from the light elements (O, Na, and Al), Fe, and the neutron-capture 
elements (Zr, La, and Eu). The stellar parameters and abundances are presented in Table 

m 



3. Results 

Even within our small sample of stars, NGC 1851 exhibits large star-to-star abundance 
variations of the light elements O, Na, and Al (see Figured]). The amplitude of the variation 
is 0.62, 1.02, and 0.58 dex for [O/Fe], [Na/Fe], and [Al/Fe] respectively. O and Na are 
anticorrelated, Na and Al are correlated, and there does not appear to be an anticorrelation 
between Mg and Al. We note that the correlations remain evident regardless of whether we 
plot [X/Fe] or log e(X). Such abundance patterns have been found in all globular clusters, 
including those with and without multiple main sequences. 

We find a mean abundance [Fe/H] = —1.27 ± 0.03 (a = 0.09) with values spanning 0.25 
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dex. The ratio [Fe/H] appears to exhibit a range that may be larger than that expected from 
observational uncertainties. It is not clear whether this is a real (though small) spread in the 
Fe abundance, or merely an artefact of small number statistics or the limited set of Fe lines 
available given the restricted wavelength coverage. Figure [T] shows that the Fe abundance is 
not correlated with the Al abundance. 

The most striking result of our analysis is that NGC 1851 has a significant star-to-star 
abundance variation of the light s-process element Zr and the heavy s-process element La. 
Figure [2] shows two pairs of stars with essentially identical stellar parameters. Therefore, 
any difference in line strengths must be due to differences in abundances between the pairs 
of stars. In this Figure, the strengths of Fe and other atomic lines are very similar between 
the pairs of stars, as expected. However, the Zr and La lines show large differences and 
therefore the Zr and La abundances must differ considerably. From a quantitative analysis, 
the amplitude of the variation is 0.56 and 0.49 dex for [Zr/Fe] and [La/Fe] respectively. This 
Figure also shows t hat star 112 is CN-strong (and La-rich) relative to star 003. Recall that 



Hesser et al.l (119821 ) found extremely CN-strong stars in this cluster and that the Sr and Ba 
lines were enhanced in these stars. Unfortunately Zr and La are the only s-process elements 
that we could measure due to the wavelength coverage and we were also unable to measure C 
or N abundances. We note that the strengths of the 6141. 7lA and 6496. 90A Ball lines differ 
from star-to-star with the Zr- and La-rich stars exhibiting stronger Ba lines. However, the 
Ba lines are saturated (> 190mA) such that reliable abundances cannot be readily measured 
and therefor e we c annot quantify the star-to-star Ba abundance dispersion. According to 



Hesser et al.l (Il982l ). star 333 is CN- normal and we find no enhancement of Zr or La. 



In Figure [3j the s-process elements Zr and La are correlated with Al and anticorrelated 
with O. Once again, the correlations are present regardless of whether we plot [X/Fe] or log 
e(X). The abundances of Zr and La cluster around two distinct values. The three La-rich 
stars have a mean value [La/Fe] = 0.61 while the five La-normal stars have a mean value 
[La/Fe] = 0.27. The three La-rich stars are also Zr-rich with a mean value [Zr/Fe] = 0.45 
and the five La-normal stars are also Zr- normal with a mean value [Zr/Fe] = 0.14. Curiously, 
one of the La-normal stars (star 209) appears to be relatively Zr-rich, [Zr/Fe] = 0.35. Setting 
aside this star, the four remaining Zr-normal stars have a mean value [Zr/Fe] = 0.09. The 
mean ratio of heavy to light s-process elements is very similar for the three La-rich, Zr-rich 
stars and for the five La-normal, Zr-normal stars, [La/Zr] = 0.16 and 0.13 respectively. 

Finally, we draw attention to the abundance of the r-process element Eu. The average 
ratio, [Eu/Fe] = +0.71, is rather hig h compared to oth er globular clusters and field stars 



at the metallicity of NGC 1851 (e.g., iPritzl et al.ll2005l ). We do not find any evidence for 



a star-to-star Eu abundance dispersion or bimodality. Eu is not correlated with the light 
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elements O or Al nor with the s-process elements Zr and La. 



Discussion 



Milone et al.l (120071 ) explored various possibilities for explaining the double subgiant 
branch in NGC 1851. If age is the sole parameter, a difference of 1 Gyr would be required. 
If composition is the sole parameter, any abundance variation invoked to account for the 
double s ubgiant branch is he avily constrained by the width of the main sequence (and giant 
branch). iMilone et al.l (120071 ) estimate that the maximum possible abundance variation on 
the main sequence would be A[Fe/H] = 0.1 dex or a helium abundance AY = 0.026. They 
argue that the He abundance alone cannot explain the observed subgiant branch. They 
suggest that a population with a 0.2 dex metallicity increase could explain the subgiant 
branch, but that such an abundance spread is precluded by the width of the main sequence 
and giant branch. However, a combination of increased [Fe/H] by 0.2 dex and helium from 
Y = 0.247 to Y = 0.30 could reproduce the observed subgiant branch behavior as well as 
ensuring a sufficiently narrow main sequence and giant branch. 

Available measurements indicate that the abundances of Fe, O, Na, Al, and neutron 
capture elements do not change as a function of evolutionary status f rom the main sequ ence 



to the giant branch in globular clusters (e.g., iGratton et al.l l200ll . iJames et al.l 120041 and 



Cohen fe Melendezll2005h . If we assume that the abundances for giant stars are representative 
of subgiant stars in NGC 1851, then our measured abundances eliminate the possibility that 
the subgiant branch split is due to an abundance difference of 0.2 dex for Fe and a shift of Y 
from 0.247 to 0.30. While our sample size is limited, there is no hint of a bimodal distribution 
of Fe abundances separated by 0.2 dex, as required for the ch emical composition explanation 
for the double subgiant branch. Indeed, IMilone et al.l (120071 ) noted that such a possibility is 
also precluded by the magnitudes of blue HB stars. Interestingly, the dispersion in our Fe 
abundances (a = 0.09 dex) is within the 0.1 dex "limit" imposed by the width of the main 
sequence. 

The abundances of additional elements add further complexity to NGC 1851. The star- 
to-star abundance variation of O, Na, and Al and the anticorrelation of O and Na exist in 
NGC 1851 and indeed in all globular clusters. While our sample size is small, the amplitude 
of these abundance variations is comparable to clusters at a similar metallicity. However, 
the light s-process element Zr and the heavy s-process element La exhibit a large star-to- 
star abundance variation. (Our spectra also indicate a large range in Ba line strengths, 
and presumably Ba abundance, but we are unable to quantify the range.) Furthermore, the 
Zr and La abundances appear to have a bimodal distribution and reinforce the idea that 
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there are two stellar populations in this cluster. (The ratio of heavy to light s-elements, 
[La/Zr] is similar for the "two populations".) La and Al (and Zr and Al) are correlated 
and so the source of the Al variation is likely to be the source of the La vari ation. The 



hint of a correlation between Al and Zr has also been found in NGC 6752 (jYong et al. 



20051 ). Asymptotic giant branch (AGB) stars may produce large enhan cements of Al and 



La, albeit via different pro cesses in AGB stars of different masses (e.g., iBusso et al.lll999 



Karakas fc Lattanziol 120031 ) . The He, C, N, and O produced by such AGB stars may play a 
role on the age differences for the subgiant branch populations. 



Milone et al.l (120071 ) raise the tantalizing prospect that the relative frequency of stars on 
the fainter /brighter subgiant branches (45% vs. 55%) roughly matches the relative frequency 
of HB stars bluer/redder than the instabil ity strip (37% vs. 63%). The fraction of CN- 
strong, Ba-strong, Sr-strong stars found by iHesser et al.l (jl982j) was ~ 40%. In this study, 
the fraction of Zr-strong, La-strong stars was also ~ 40%. Based on the relative numbers, 
we speculate that the brighter subgiant branch stars have "normal" CN, Sr, Zr, Ba, and La 
abundances and populate the red HB. We also speculate that the fainter subgiant branch 
stars are CN-strong, enriched in Sr, Zr, Ba, and La, and populate the blue HB. Indeed if 
the CN, Zr, and La enriched stars are also He-rich, they are expected to populate the blue 
HB f 1Sweigartlll997l ). However, we caution that the relative numbers of the two populations 
will not be the same at all evolutionary phases. Nevertheless, the abundances in subgiant 
branch, giant, and red HB stars can test this hypothesis. 

Since the fainter subgiant branch stars are older than the brighter subgiant branch 
stars (assuming age is the only parameter controlling the double subgiant branch), then our 
tentative speculation would require that the Zr-rich, La-rich stars are younger than the Zr- 
normal, La-normal st ars. While this unpalatable scenario would be a challenge to explain, 
Villanova et al.l (120071 ) find that the most metal-rich subgiant stars in u Centauri are among 
the oldest stars in the cluster. Clearly, the sequence of events that led to the formation of 
the different stellar populations in NGC 1851 and uj Centauri requires some imagination. 
On the other hand, an analysis of a larger sample of giant stars may reveal that the CN- 
strong, Zr-rich, La-rich stars are associated with the majority population of brighter subgiant 
branch stars. In this case, if the CN-strong, Zr-rich, La-rich stars are the progeny of the 
CN-normal, Zr-normal, La-normal stars, the assumed 1 Gyr age difference would then place 
constraints upon the mass range of the AGB stars that produced the CN, Zr, and La excess. 
At the metallicity of NGC 1851, Karakas (2007 priv. comm.) suggests that the minimum 
mass (ZAMS) of the AGB stars would be approximately 2 M . Further, the AGB stars that 
produce La have M > 1 M Q and the AGB stars that produce Al have M > 4 M Q . 



Finally, we note that the only other globular cluster that exhibits a large range in C, N, 
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O, Na, Al, Sr, Zr, Ba, and La abundances is u Centauri. Both NGC 1851 and ui Centauri 
display multiple subgiant branches which are presumably due to stellar populations with 
different ages and compositions. Further, NGC 1851 and ui Centauri are the only clusters 
that display a large variation in the Stromgren mi index, traditionally used as a metallicity 
indicator. NGC 1851 appears to exhibit a bimodal m,\ distribution on the giant branch 
suggesting an extreme variation in C and/or N abundances (Grundahl 2007 in preparation). 
In the context of the formation of the Galactic globular cluster system, we speculate that 
NGC 1851 may represent a "bridge" between to Centauri (large variation of all elements) 
and NGC 6752-like clusters (constant Fe but large variations of light elements C-Al). A 
more detailed analysis of additional elements and additional stars in NGC 1851 is of great 
interest. 
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supported by the Danish National Science Research Council. This research was supported 
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Table 1. Stellar parameters and abundances 



Star a 


T eff (K) 


log g (cgs) 


Vt (km s l ) 


[Fe/H] 


[O/Fe] 


[Na/Fe] 


[Mg/Fe] 


[Al/Fe] 


[Zr/Fe] 


[La/Fe] 


[Eu/Fe] 


003 


4075 


0.60 


1.80 


-1.28 


0.52 


0.13 


0.25 


0.32 


0.02 


0.28 


0.83 


095 


4025 


0.15 


1.95 


-1.40 


0.44 


0.23 


0.50 


0.46 


0.14 


0.30 


0.75 


112 


4025 


0.25 


1.85 


-1.34 


0.13 


0.74 


0.42 


0.76 


0.38 


0.69 


0.69 


151 


4175 


0.65 


1.80 


-1.32 


0.61 


-0.14 


0.35 


0.26 


0.09 


0.27 


0.82 


209 


4400 


1.15 


1.75 


-1.20 


0.49 


-0.27 


0.40 


0.21 


0.35 


0.30 


0.70 


329 


3875 


0.10 


1.55 


-1.20 


0.39 


0.34 


0.33 


0.49 


0.58 


0.60 


0.65 


333 


4225 


0.75 


1.70 


-1.25 


0.44 


-0.21 


0.36 


0.18 


0.10 


0.20 


0.65 


395 


3975 


0.45 


1.65 


-1.15 


-0.01 


0.62 


0.39 


0.63 


0.40 


0.55 


0.60 



l Star names taken from IStetsonl (jl98ll ). 
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Fig. I.— Clockwise from upper left, [Na/Fe] vs. [O/Fe], [Mg/Fe] vs. [Al/Fe], [Fe/H] vs. 
[Al/Fe], and [Al/Fe] vs. [Na/Fe]. Linear least squares fits to the data are shown (slope and 
associated error are included). O and Na are anticorrelated and Al and Na are correlated. 
Mg and Al do not appear to be correlated nor are Fe and Al. 
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Fig. 2. — The upper panel shows stars 095 and 395 centered on the Zrll lines. The stellar 
parameters are very similar as confirmed by the atomic lines of similar strength. However, 
the Zr lines differ considerably. The lower panel shows stars 003 and 112 centered on the La II 
line. The stellar parameters are very similar as are the strengths of atomic lines. However, 
the La line differs considerably as do the CN molecular lines. 
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Fig. 3. — Clockwise from upper left, [Zr/Fe] vs. [Al/Fe], [La/Fe] vs. [Al/Fe], [La/Fe] vs. 
[O/Fe], and [Zr/Fe] vs. [O/Fe]. Linear least squares fits to the data are shown (slope and 
associated error are included). Zr and La increase with increasing Al and decrease with 
increasing O. 



